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Resistance Noise Near to Electrical Breakdown: 
Steady State of Random Networks as a Function of the Bias 
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A short review is presented of a recently developed computational approach which allows the 
study of the resistance noise over the full range of bias values, from the linear regime up to electrical 
breakdown. Resistance noise is described in terms of two competing processes in a random resistor 
network. The two processes are thermally activated and driven by an electrical bias. In the linear 
regime, a scaling relation has been found between the relative variance of resistance fluctuations and 
the average resistance. The value of the critical exponent is significantly higher than that associated 
with 1/f noise. In the nonlinear regime, occurring when the bias overcomes the threshold value, the 
relative variance of resistance fluctuations scales with the bias. Two regions can be identified in this 
regime: a moderate bias region and a pre-breakdown one. In the first region, the scaling exponent 
has been found independent of the values of the model parameters and of the bias conditions. A 
strong nonlinearity emerges in the pre-breakdown region which is also characterized by non-Gaussian 
noise. The results compare well with measurements of electrical breakdown in composites and with 
electromigration experiments in metallic lines. 
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I. INTRODUCTION 

A large interest exists in the recent literature concern- 
ing the electrical and mechanical stability of different 
physical systems^^^^. In fact, from a technological point 
of view, the increasing level of miniaturization of elec- 
tronic devices enhances the importance of degradation 
and failure processes^, ^'',^^~^^. On the other hand, from 
a theoretical point of view, the study of the response 
of disordered systems to high external stresses, implied 
in the study of breakdown and fracture processes, is of 
great help to understand the properties of these sys- 
tems and, more generally, to develop the knowledge of 
nonequilibrium systems^"^^,"^^""^^. It is well known that 
the application of a finite stress (electrical or mechani- 
cal) to a disordered material generally gives a nonlinear 
response, which ultimately leads to an irreversible break- 
down (catastrophic behavior) in the high stress limit ^~^. 
Such breakdown phenomena have been successfully stud- 
ied by using percolation theories^~^,^~^,^^~^^,^°,^^~''^. In 
particular, by focusing on electrical breakdown, a large 
attention has been devoted to the determination, by 
both theory^"3,^''^'2°,3^"^2 a,nd experiments^"^, ^^"^^ of 
the critical exponents describing the resistance and the 
variance of resistance fluctuations in terms of the con- 
ducting particle (or defect) concentration"^^'"^^. In fact, 
it is well known that the study of the resistance fluctua- 



tions is a fundamental tool to extract information about 
the system stability ^~'*,^^~^^,'^^~^*. In spite of the wide 
literature on the subject few attempts have been made so 
£g^j.2i,22 -j-g describe the behavior of a disordered medium 
over the full range of the applied stress, i.e. by studying 
the response of the system to an external bias when the 
bias strength covers the full range of linear and nonlin- 
ear regimes. On the other hand, important information 
are expected from this kind of study, like: precursor phe- 
nomena, role of the disorder, existence of scaling laws, 
predictability of breakdown, etc. 

Recently, a new approach has been proposed to in- 
vestigate the resistance noise in thin films over the full 
range of bias values, from the linear regime up to electri- 
cal breakdown**^"^^. In this paper, I will present a brief 
review of this approach. The model consists of describ- 
ing the resistance noise in a conducting film in terms of 
two competing processes taking place in a random resis- 
tor network (RRN)"'^"^^. The two processes are driven 
by the joint effect of the electrical bias and of the heat 
exchange with a thermal bath. The breaking of a single 
resistor of the network is taken as a simple and general 
manner to account for different microscopic mechanisms 
which are responsible for the degradation of the electri- 
cal properties of a small region of the film. This defect 
generation process is taken to occur in competition with 
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an opposite process, named defect recovery, which mim- 
ics the healing mechanisms. Electromigration of metaUic 
lines'*'^^ instabihty of the electrical properties of compos- 
ites or semicontinuous metal films^'^,^'^"^^,^^'^'^ or soft di- 
electric breakdown of ultra-thin oxides^'^^, are examples 
of phenomena that can be successfully described by this 
approach^^''^^'^^"^^. In the case of electromigration phe- 
nomena, for example, the defect generation corresponds 
to the formation of voids induced by the electronic wind, 
while the defect recovery is related to the void healing 
due to mechanical stress and thermal gradients inside a 
metallic film'''^^. In the case of composites, the electri- 
cal instability arises from the opening of new conductive 
channels as a consequence of thermally activated hopping 
processes^'^^. On the other hand, the same processes can 
also lead to a suppression of already existing channels. 
In the case of soft dielectric breakdown of ultra-thin ox- 
ides it is instead the competition between electron trap- 
ping and detrapping which plays a key role on leakage 
currents'''^^''^''. 

All these phenomena and many others^""^,^^'^^ can be 
modeled in terms of competition between two opposed 
processes taking place in a RRN. Monte Carlo simula- 
tions show that, depending on the bias strength, an irre- 
versible failure or a stationary state of the RRN can be 
achieved. By focusing on the steady-state, the behavior 
of the average network resistance and the properties of 
the resistance fluctuations are analyzed as a function of 
the bias. At low bias, an effective defect generation prob- 
ability can be defined controlling the network behavior'*^. 
In this Ohmic regime, a scaling relation has been found 
between the relative variance of resistance fluctuations 
and the average network resistancc^^. The properties 
of the nonlinear regime, occurring when the bias over- 
comes a threshold value, are studied for different values 
of the model parameters and for different bias conditions 
(constant voltage or constant current)^"'^^. Scaling re- 
lations are found relating the average resistance and the 
resistance noise with the external bias. In particular, 
a strong nonlinearity emerges in the prc-breakdown re- 
gion which is also characterized by a non-Gaussian noise. 
The results compare well with measurements of electri- 
cal breakdown in composites^'^^'^^ and in semicontinuous 
metal films^^ and with electromigration experiments in 
metalhc lines^^'^^. 



II. THEORY 

The RRN consists of a two-dimensional square-lattice 
network made of Ntot resistors of resistance r^'^^ . Differ- 
ent geometries can be studied; here we consider the sim- 
plest one which is the square N x N, where N determines 
the linear size of the network. In practical applications 
the value of N can be related to the ratio between the size 
of the film and that of the characteristic grain. The RRN 
exchanges heat with a thermal bath at temperature Tg 



(substrate temperature). A current /, kept constant, is 
applied to the network through electrical contacts at the 
left and right hand sides. Alternatively, constant volt- 
age conditions can be considered. The current flowing in 
the n-th resistor of the network is denoted by i„. The 
resistances are taken linearly dependent on temperature 
through a temperature coefficient of resistance, a: 



rn{T„) = ro[l -I- a(T„ - Tq)] 



(1) 



The local temperature r„ is calculated by adopting the 
biased percolation model^^'^^ as: 
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where Nndg is the number of first neighbors around the 
nth resistor. The parameter A, measured in (K/W), de- 
scribes the heat coupling of each resistor with the thermal 
bath and it determines the importance of Joule heating 
effects. The parameter B is taken to be equal to 3/4 to 
provide a uniform heating in the perfect network configu- 
ration. It must be noticed that Eq. (2) implies an instan- 
taneous thcrmalization of each resistor at the value r„, 
thus, by adopting Eq. (2), we are neglecting for simplicity 
time dependent effects which are discussed in Ref.^. 

Here, the initial state of the network is chosen coin- 
ciding with the perfect network (r„ =7^0) of resistance 
Rq. Nevertheless, a disordered initial configuration can 
easily be introduced in the present approach to analyze 
the role of the initial disorder on the RRN response. 
We assume that two competing processes act to deter- 
mine the RRN evolution. The first process consists of 
generating fully insulating defects (resistors with very 
high resistance, i.e. broken resistors) with probability^*^ 
Wo.n — Gxp{~Eo/KBTn), whcrc Eo is a characteristic 
activation energy and Kb the Boltzmann constant. The 
second process consists of recovering the insulating de- 
fects with probability Wfl,„ = exp{—Eii/KBTn), where 
Eji is an activation energy characteristic of this second 
process^^'"*^. Thus, the first process consists in a percola- 
tion of broken resistors within a network of active resis- 
tors. This percolative process is contrasted by a recovery 
process. This second process can also be seen as a per- 
colation of active resistors within an insulating network. 
For A ^ 0, Eq. (2) implies that both the processes (de- 
fect generation and defect recovery) are correlated pro- 
cesses. Indeed, the probability of breaking (recovering) 
a resistor is higher in the so called "hot spots" of the 
RRN^e. On the other hand, for ^ == Eq. (2) yields 
T„ = To, which corresponds to the competition of two 
random processes^^'^^. The same is true for vanishing 
small bias values, when Joule heating effects are negligi- 
ble. 




FIG. 1. Evolution of 40 x 40 network obtained for a value 
of the external current / — 1.2 (A). The different colors, from 
black to violet, correspond to increasing values of Vn from 1 
to 15 (Q). The pattern shown corresponds to the iteration 
step t = 800. 
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FIG. 3. Resistance evolutions for increasing bias values. 
Going from bottom to top, the curves are obtained for 
/ = 0.01, 0.1, 0.5, 1.0, 1.5, 1.8, 2.1 (A) and they correspond to 
RRN steady states. Precisely, the highest curve shows an un- 
steady evolution of the resistance obtained for I = 2.2 (A), 
the curve just below corresponds to the breakdown current 
h = 2.1. 




FIG. 2. Pattern corresponding to the same evolution of 
Fig. 1 and to the iteration step t = 805. 



The sequence, defect-creation and defect-recovery, is 
then iterated. Depending on the material parameters 
{Eo, Efi, A, a, tq, N) and on the external conditions 
(specified by the kind of applied bias and by the bath 
temperature) a steady state or an irreversible breakdown 
characterized by a critical fraction of defects Pc (per- 
colation threshold) are reached. In the first case, the 
network resistance fluctuates around an average value 
< R >. In the second case, R diverges due to the exis- 
tence of at least one continuous path of defects between 
the upper and lower sides of the network"^^. We note 
that in the limit of a vanishing bias (two random pro- 
cesses) and infinite lattices (N -^ oo), the expression: 
Er < Ed + KbTo ln[l + eyi^{- E d / K bTq)] provides a 
sufficient condition for the existence of a steady state*^. 
In the general case of biased processes, no analytical ex- 
pressions are available for the steady state condition and 
thus it is necessary to rely on numerical simulations only. 



The evolution of the RRN is obtained by MC simu- 
lations carried out according to the following procedure, 
(i) Starting from the perfect lattice with given local cur- 
rents, the local temperatures r„ are calculated according 
to Eq. (2); (ii) the defects are generated with probabil- 
ity Wd and the resistances of the unbroken resistors are 
changed as specified by Eq. (1); (iii) the currents «„ are 
calculated by solving Kirchhoff's loop equations by the 
Gauss elimination method and the local temperatures are 
updated; (iv) the defects are recovered with probability 
Wr and the temperature dependence of unbroken resis- 
tors is again accounted for; (v) R, i„ and r„ are finally 
calculated and the procedure is iterated from (ii) until 
one of the two following possibilities is achieved. In the 
first, the percolation threshold is reached. In the second, 
the RRN attains a steady state; in this case the iteration 
runs long enough to allow a fluctuation analysis to be 
carried out. Each iteration step can be associated with 
an elementary time step on an appropriate time scale (to 
be calibrated with experiments). In this manner it is 
possible to represent the simulated resistance evolution 
over either a time or a frequency domain according to 
convenience. Except when differently specified, the sim- 
ulations are performed by taking the following values for 
the parameters, which are chosen as reasonable values: 
N = 75, ro ranges from 1 ^ 10 (f7), a = IQ'^ (K'^), 
A = 5 X 10^ (K/W), Ed = 0.170 (eV), Er in the range 
0.026 H- 0.155 (eV) and To = 300 (K). The values of the 
external bias range from 0.001 < I < 3.0 (A) under con- 
stant current conditions, and from 0.001 < V < 3.5 (V) 
under constant voltage conditions. Further details about 
numerical simulations can be found in Refs.^'*''*^"^^. 
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FIG. 4. Average resistance versus the external current. 
The resistance is normalized to the linear regime value. 
The different curves are obtained for Ed = 0.17 (eV) and 
Eh, = 0.026 (eV), while ro ranges from 1 to 10 (fi). 



states of the RRN, the upper curve to an evolution to- 
ward electrical breakdown. In particular, it must be no- 
ticed that the first two curves (black and cyan) , obtained 
for / = 0.01 (A) and / = 0.1 (A), are practically overlap- 
ping. By contrast, all the other steady state curves fluc- 
tuate around signiflcantly different average values with 
increasing fluctuation amplitudes. Thus, the figure well 
illustrates a general feature of the model: the ability of 
describing both the linear and the nonlinear regimes up 
to the breakdown. 
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FIG. 5. Average resistance versus the external current. 

The resistance is normalized to the linear regime value. The 

different curves are obtained for Ed ~ 0.17 (eV) and ro = 1 

(n), while Er ranges from 0.026 to 0.155 (eV). 
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FIG. 6. Normalized average resistance versus the external 
bias: constant current for the green and blue curves, constant 
voltage for the red and black ones. Precisely, the green and 
red curves are obtained for a = 10~ (K~ ) while the blue 
and black for a — 0. The other parameters are specified in the 
text. The bias values are in A and in V units under constant 
current and constant voltage, respectively. 



III. RESULTS 



Figures 1 and 2 display the evolution of a 40 x 40 net- 
work, made of resistors with rg = 1.0 (fi), characterized 
by a value of En = 0.043 (eV) and subjected to a bias 
current / = 1.2 (A). For this choice of the parameters, 
the network becomes unstable for currents greater than 
lb = 1.1 (A). Therefore, in the case considered, the net- 
work evolves towards the electrical breakdown which is 
reached after 809 iteration steps. Precisely, Fig. 1 shows 
the formation of filamented clusters of defects (missing 
resistors) characteristic of biased percolation^'*'^^. Figure 
2 reports the pattern very near to the final breakdown. 
The incipient cluster of defects, perpendicular to the di- 
rection of the applied current is well evident, together 
with the formation of hot spots. This kind of damage pat- 
tern perfectly agrees with the pattern observed in metal- 
lic lines failed as a consequence of electromigration*'^^. 

Typical resistance evolutions of a 75 x 75 RRN, with 
the same parameters used for Figs. 1 and 2, are reported 
in Fig. 3 for different values of the bias current. The 
first six curves from the bottom correspond to steady 
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FIG. 7. Log-log plot of the relative variation of the average 
resistance versus the normalized current. Only the nonlinear 
regime is shown. The data correspond to different values of 
Er and to different values of ro and they are the same shown 
in Figs. 4 and 5. The dashed line fits the data in the moderate 
bias region. 
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Log-log plot of the relative variation of the average 
resistance versus the normalized current. Only the nonlinear 
regime is shown. The data indicated with green stars are ob- 
tained for a = 10~^ (K~^), while the data indicated with blue 
squares for a = 0. The dashed black lines fit the data with a 
power-law of exponent 2.1, the magenta line with a power-law 
of exponent 3.1 and the black solid line corresponds to an ex- 
ponent 3.7. 
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FIG. 9. Log-log plot of the relative variation of the average 
resistance versus the normalized voltage. The data indicated 
with red triangles are obtained for a = 10~^ (K~^), while 
the data indicated with black diamonds for a = G. The black 
lines fit the data with a power-law of exponent 2.1 while the 
magenta line with a power-law of exponent 2.8. 



Details about the criteria used for the determinations of 
/o and lb can be found in Ref.^^. It must be noticed that 
by increasing the initial RRN resistance, both /;, and /q 
decrease but the ratios /b/Io and < R >i, / < R >q are 
found to be independent of the initial resistance value^°. 
Moreover, these ratios are found independent of the net- 
work sizes^'^. Both these behaviors agree with electrical 
breakdown measurements performed in composites^^ and 
in semicontinuous metal films^^. 



Figure 5 reports the normalized average resistance as 
a function of the bias current for different values of the 
recovery activation energy, specified in the figure. All the 
curves are now obtained by taking rg = 1 (fJ). In this 
case, not only Ii, and Iq decrease with Efi but also the 
ratios h/Io and < R >b / < R >o diminish^". Thus, 
in contrast to what happens for the effect of the initial 
resistance, the decreasing robustness of the system is as- 
sociated with a reduction of the extent of the nonlinear 
region. It must be underlined that if the value of En 
is sufficiently different from its maximum value Eh^max 
(determined by the stability condition in the vanishing 
bias Hmit, discussed in Sec. 2) the threshold current /& 
is associated with a first order transition. In fact, the 
"last" stable state of the RRN, corresponding to /{,, is 
characterized by an average defect fraction < p >{, which 
is smaller than pc- Therefore, the correlation length and 
the correlation time remain finite at the threshold^^ , de- 
noting a non critical behavior. It must be noticed that 
the prediction of a first order transition is confirmed by 
experiments^'^^. 



Figure 4 shows the average values of the RRN resis- 
tance as a function of the bias current. Each value is 
calculated by considering the time average on a single 
steady-state realization and then averaging over 20 in- 
dependent realizations. Further details can be found in 
Ref.^"'^. The different curves correspond to different val- 
ues of ro , as specified in the figure, and to a common value 
of En = 0.026 (eV). At the lowest biases the average re- 
sistance is independent of the bias values and it takes a 
value < R >o, which represents the linear response of the 
network (Ohmic regime)'*^. For this reason in Fig. 4 the 
average resistance has been normalized to < i? >o. At 
increasing biases, when the current overcomes a certain 
value, /q, the average resistance starts to become depen- 
dent on the bias. Thus, Iq sets the current scale value 
for the onset of nonlinearity. The average resistance in- 
creases with bias up to a value < i? >{,, which corre- 
sponds to a threshold current h- Above this threshold 
current the RRN undergoes an irreversible breakdown. 
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FIG. 10. Percolation threshold versus the external current 
for difi'erent values of the recovery energy Er which ranges 
from 0.026 (eV) to 1.0 (eV). 
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FIG. 11. Relative variance of resistance fluctuations ver- 
sus the external current. Magenta, green and cyan curves 
correspond to decreasing values of Er. The black curves are 
obtained for the same value of Er used for the cyan curve and 
they differ for the value of ro. Precisely, cyan circles ro — 1.0, 
open square ro = 2.5, black triangles ro = 5.0, open triangles 
ro = 7.5, crosses ro — 10.0. The values of ro are expressed 
in Q. The red curve with down triangles is obtained for the 
same parameters of the cyan curve but for a smaller network, 
of size A*' = 36. 
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FIG. 12. Relative variance of resistance fluctuations versus 
I < p > — Pel in the linear regime. This regime corresponds 
to the competition of two random processes. 

The effect of the temperature coefficient of the resis- 
tance, a, is shown in Fig. 6, where we have also reported 
the behaviors of the average resistance calculated under 
constant current and under constant voltage conditions. 
The simulations are performed by taking ro = 1 and 
En = 0.043 (eV). Similarly to the constant current case, 
also under constant voltage two threshold voltage values 
exist: Vq and Vf,, corresponding to the nonlinearity onset 
and to the electrical breakdown, respectively. The ratio 
Vb/Vo is found to be significantly higher than the ratio 
h/Io^^- This fact reflects the greatest stability of the 
system when biased under constant voltage than under 
constant current. This property is further emphasized 
by the fact that the increase of the resistance in the pre- 
breakdown region exhibits a lower slope under constant 
voltage than under constant current conditions, as dis- 
cussed below. It must be noticed that in spite of the 
significant difference of the ratios Vb/Vo and h/Io, the 
ratio < R >b / < R >Q remains practically the same 
under the different bias conditions'^, in agreement with 
experiments^^. 
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FIG. 13. Log-log plot of the relative variation of E versus 
I/Io- The data shown are obtained for different values of Er, 
ro and A^ and they are the same of Fig. 11. The black line 
fits the different sets of data in the moderate bias region with 
a power-law of exponent 2.0, while the blue line of slope 5.0 
in the high bias region is only for a guide to the eyes. 

To analyze the dependence on the bias of the aver- 
age resistance and to investigate the existence of scaling 
relations and their universality'^^. Fig. 7 reports on a log- 
log plot the relative variation of the average resistance, 
A < R> / < R >o= {< R> - < R >o)/ < R >o, as 
a function of I/Iq for different values of tq and Ej^. The 
data are the same of Figs. 4 and 5. The figure shows 
that all the curves collapse onto a single one and that 
the relative variation of < i? > scales with the ratio I/Iq 
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= g{I/Io), 9(1/ lo) ^ 1 + a(I/Iof (3) 



with the exponent 6 = 2.1 ± 0.1 and a a dimension- 
less coefficient. However, it must be noticed that a su- 
perquadratic behavior of < i? > emerges in the pre- 
breakdown region, where the relative variation of < i? > 
is characterized by a power law (I/Iq)^' with an expo- 
nent 6 J = 3.7 ± 0.3^^. Furthermore, this superquadratic 
behavior emerges only for RRN sufficiently robust. In 
fact, at increasing the recovery activation energy Er, the 
ratio Ib/Io becomes smaller, i.e. the stability region is re- 
duced and the dependence on the current of the average 
resistance can remain quadratic over the entire region. 
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FIG. 14. Relative variance ol resistance fluctuations versus 
the external bias: constant current for the green and blue 
curves, constant voltage for the red and black ones. Precisely, 
the green and red curves are obtained for a = 10~3 (K~^) 
while the blue and black for a = 0. The other parameters are 
specified in the text. The bias values are in A and in V units 
under constant current and constant voltage, respectively 
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FIG. 15. Log-log plot of the relative variance of resistance 
fluctuations versus the normalized current. The data indi- 
cated with green stars are obtained for a = 10~3 (K~^), 
while the data indicated with blue squares for a = 0. The 
black dashed lines fit the data with a power-law of exponent 
2.1, the black solid line with a power- law of exponent 5.4, 
while the magenta solid line corresponds to an exponent 4.0. 
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FIG. 16. Log-log plot of the relative variance of resistance 
fiuctuations versus the normalized voltage. The data indi- 
cated with red triangles are obtained for a = 10~3 (K~^), 
while the data indicated with black diamonds for a = 0. The 
black dashed lines fit the data with a power-law of exponent 
2.0, the black solid line with a power-law of exponent 3.0, 
while the magenta solid line corresponds to an exponent 4.0. 





4 


r 


\ 1=1 -8 (A) 




2 


i 


\o 
























f 


\ottP ' 


r 






\0 ° 












-2 




\ o 




-4 
p. 




CD 



-0.02 



0.02 



FIG. 



R-<R> (ii) 
17. Distribution function of the resistance fluctua- 



tions for an applied current of I 

rameters are specified in the text 

therefore the dashed curve corresponds to a Gaussian distrl 

bution. 



= 1.8 (A), the other pa- 
The scale is a linear-log, 



lO'' 
10' 
10° 

n 

E10 
10-^ 
10"' 
10" 




10 



100 



1000 



f(a.u) 



FIG. 18. Power spectral density of resistance fiuctuations 
under constant current conditions. The solid curve is obtained 
for / = 1.5 (A), the dotted one for / = 1.8 (A). 

In the log-log plot reported in Fig. 8 we compare the 
relative variations of the resistance versus the ratio I /Iq 
calculated for two different values of the temperature 
coefficient: a = 10"'^ (K~^) and a = (the data are 
the same as those shown in Fig. 6 with stars and open 
squares respectively). We can see that the two sets of 
data do not collapse onto the same curve, i.e. they do not 
belong to the same universality class. Moreover, the ef- 
fect of a nonzero value of a is to significantly enhance the 
superquadratic dependence of the resistance on the bias 
under constant current conditions. By contrast, the op- 
posite is true under constant voltage conditions, as shown 
in Fig. 9, where the data obtained for a — 10"'^ (K~^) 
are again compared with those obtained for a = (same 
data shown in Fig. 6 with triangles and full squares, 
respectively). In fact, in this case, the effect of a is to 
depress the superquadratic behavior. Remarkably, under 
constant voltage conditions, the average resistance scales 
quadratically with the applied voltage over the whole re- 
gion of bias values up to breakdown, as^^: 



<R> 



V 



<R>o 



= 1 



(4) 



where a' is a dimensionless coefficient and the value of 9 
is the same of that in Eq. (3). 



The quadratic dependence on the bias in the moderate 
bias region, common to aU the data shown in Figs. 7, 8 
and 9, must be considered a general feature of the model, 
in agreements with experiments^^ and can be explained 
as follows. It is well known^^ that the resistance of suf- 
ficiently large RRNs subjected to random percolation is 
related to the fraction of broken resistors by the scaling 
relation: 



R^\p-Pc 



(5) 



where Pc = Pco — 0.5 for a square lattice and /x takes a 
universal value, known from very accurate calculations: 
fj, = 1.303 for two-dimensional RRNs'^^. This scaling 
relation should hold also for steady states of RRNs re- 
sulting from the superposition of two opposite random 
processes. In this case, Eq. (5) relates the average net- 
work resistance < R > with the average fraction of de- 
fects < p >. Numerical simulations reported in Ref. 
confirm this statement. On the other hand, when a sin- 
gle biased percolation is present Pc becomes a function 
of the bias strength^^'^^ (see black circles in Fig. 10), 
moreover the value of /i is no more universal and it de- 
pends on the biasing conditions''*'^'*. When two compet- 
ing biased processes are present, leading to steady states 
of the RRN, the use of Eq. (5) becomes quite problem- 
atic. In fact, in this case pc depends also on the re- 
covery activation energy, as shown in Fig. 10. In this 
figure the black circles represent Pc versus the bias when 
En » Ed so that no recovery of defects occurs. In this 
case the RRN is unstable for any value of the applied cur- 
rent and, at vanishing bias, Pc takes the random percola- 
tion value Pco- At increasing bias, the correlated growth 
of defects along filamented clusters reduces pc- At still 
higher values of the bias, a multi-channel filamentation 
emerges'^ which makes pc an increasing function of /. 
When Eli < Ed a stability region appears for / < /t, 
whose extent increases at decreasing Eji, as shown in 
Fig. 5. Therefore, the corresponding curves in Fig. 10 
start at progressively higher values of /. In spite of these 
complications, at moderate bias and when < p ><C Pc, 
by truncating Eq. (5) to the first order in < p > /pc, we 
can write: A < R >« C{fio,Pco) A < p >. Further- 
more, in the spirit of a mean-field theory, it is possible 
to see'^ that the relative variation of the average defect 
fraction: A<p>/<p>ooc (I/Iq)^. Thus, also 
A < i? > oc {I /IqY . At high bias, terms of order higher 
than the first in the expansion of < i? > in terms of 
< p > Ipc and the dependence on the bias of Pc become 
important. Furthermore, A < p > is no more a quadratic 
function of the bias'^. Therefore, the average resistance 
acquires the superquadratic behavior shown in Figs. 7, 
8 and 9. On the basis of the previous results and of the 
above discussion, the behavior of the average resistance in 
the moderate bias region is expected to be a universal fea- 
ture. Of course a final validation of this statement would 
require a study of the effect of the network topology (non 
square lattices, etc.) and of the boundary conditions. 



The relative variance of resistance fluctuations, S =< 
A.R? > I < R >^, is reported in Fig. 11 as a function of/ 
and for different tq and Eb. values. The data correspond 
to the same simulations shown in Figs. 4 and 5 and they 
are obtained by using the same procedure of time averag- 
ing over a single simulation and then ensemble averaging 
over 20 realizations. Only the three lowest values of E^, 
used in Fig. 5 have been reported in Fig. 11 and with 
the same symbols. In addition, the red curve with down 
triangles shows the relative variance S obtained for a net- 
work of size 36 x 36. All the other parameters are the 
same as those ones used for the curve with cyan circles 
{Eji ~ 0.026 eV and ro = 1.1^). At low bias, the relative 
variance of resistance fluctuations is found to achieve a 
constant value Sq which represents an intrinsic property 
of the system. On the other hand, a strong increase of S 
characterizes the nonlinear regime occurring for I > Iq. 

By focusing on Sq, we can see in Fig. 11 that this quan- 
tity is independent of ro but strongly dependent on the 
recovery activation energy E^. In fact Sq, corresponding 
to the vanishing bias limit, arises from the competition 
of two random processes. The study of the resistance 
noise in RRNs subjected to two random processes, has 
been performed in Ref."*^ and it pointed out the follow- 
ing results. So exhibits two regimes as a function of 
I < p > — Pco I (see Fig. 12). A nearly perfect network 
regime occurs when < AR^ >o is less than Rq/2N'^, 
while a disordered network regime appears in the oppo- 
site case. In the first case, the resistance noise is directly 
proportional to the fraction of defects'*^. In the other 
case (disordered network), the breaking and the recover- 
ing of the backbone resistors"^^ result in an enhancement 
of the resistance noise. It is noteworthy that the first 
regime is a finite size effect which disappears in the limit 
of an infinitely large network*^. In the second regime the 
data follow closely the scaling relation: 



<p> - PcO 



i-fe 



(6) 



with k = 3.1*^. A similar scaling relation, with a dif- 
ferent critical exponent kf = 1.12, has been found be- 
tween the relative variance of resistance fluctuations and 
I < p > — Pco\ in the case of 1/f noise in random 
networks"^^. On this respect, it must be noticed that 
the expression 1/f noise is used as a shorthand for a 
noise (i) spatially uncorrelated, (ii) statistically uniform 
and (iii) sufficiently small on each elementary resistors 
of the network to be treated only to the lowest order. 
Thus, it is this last feature (iii) which makes the cru- 
cial difference between the noise in Ref.^^ and the Sq 
noise given by Eq. (6), which instead arises from the 
random switch-off and switch-on of the elementary re- 
sistors. As a consequence, the noise in Eq. (6) is strongly 
sensitive to the average defect fraction. By combining 

So - I < J5 > - PcO r'' with < R >o- \<p> -PcO I"'', 
one obtains in the disordered network regime: 



^<R>f 



(7) 



with s = k/jj, — 2.6. A more detailed analysis of the 
properties of Sg can be found in Ref.^^. 

Now, by considering the behavior of S in the nonlinear 
regime, it is convenient to report its relative variation, 
(S — So)/So, on a log-log plot as a function on I/Iq. 
This is done in Fig. 13, where the same data of Fig. 11 
are shown. Similarly to what happens for the average 
resistance, two regions can be identified in the nonlinear 
regime. At moderate bias, S scales quadratically with 
the current, as: 



^0 



/(///o)^l + c(///o)'' 



(8) 



with T] = 2.0 ± 0.1 and c a dimensionless coefficient. At 
higher biases, in the pre-breakdown region, a strong su- 
perquadratic dependence appears^"'^'^. It must be re- 
marked that the behavior of S in the pre-breakdown re- 
gion is found to be dependent on Ej^. In fact, in this 



region: 



So 



j\Vl 

To, 



(9) 



with 4.5 < rji < 5.5, depending on En. Therefore, the 
straight-line with slope 5 in Fig. 13 has been drawn for 
a guide to the eyes. This behavior of the relative vari- 
ance of resistance fluctuations in the pre-breakdown re- 
gion is different from that obtained in the same region 
for the average resistance. In fact, as already discussed 
in connection with Fig. 7, Oi is found to be practically 
independent of the Eji value. The strong nonlinearity 
and nonuniversality of S in the pre-breakdown region 
has been observed in electrical noise measurements in 
conducting polymers^^ and in other materials^'^^'"*^'^*. 

The effect on S of a and of the different bias condi- 
tions (constant current or constant voltage) is considered 
in Fig. 14. The values of S in these figures are obtained 
from simulations performed for En — 0.043 (eV) and 
To = 1 (fl) (as in Figs. 6, 8 and 9). The data indicated 
by stars and open squares are obtained under constant 
current conditions with a = 10""^ (K^^) and a = re- 
spectively, while data denoted with triangles and black 
diamonds are obtained under constant voltage conditions 
for the two different values of a. Fig. 14 shows that the 
linear regime value, Sq, is independent of a and bias 
conditions. By contrast, a significantly different behav- 
ior of S is found in the nonlinear regime. In particular, 
the increase of S in the pre-breakdown region is steeper 
under constant current than under constant voltage con- 
ditions. This feature is ascribed to the better stability of 
the RRN under constant voltage, since resistance fluctu- 
ations in excess over the mean value are damped in this 
condition while they are enhanced under constant current 
conditions. The log-log plots of the relative variation of 
S as a function of I/Iq and V/Vq, reported in Figs. 15 
and 16 respectively, display the following features. The 
quadratic dependence of E on the bias in the moderate 
bias region is a feature independent of the values of a and 



of the bias conditions. Further investigations, concerning 
the effect of the network structure and of the geometry of 
the electrical contacts on this feature are necessary to es- 
tablish the universality of this behavior. By contrast, the 
superquadratic behavior in the pre-breakdown region is 
strongly influenced by both the a value and the bias con- 
ditions. Precisely, the simulations in Fig. 15, performed 
under constant current conditions with a = and with 
a — 10^"^ (K^^), show a strong enhancement of the rel- 
ative variance of resistance fluctuations induced at high 
bias by a positive value of the temperature coefficient of 
the resistance. On the other hand, by comparing the su- 
perquadratic behavior obtained, with a = 10~^ (K~^), 
under constant current conditions and under constant 
voltage, we can see that the value of rji is drastically 
reduced from 5.4 to 3.0 under constant voltage. An ex- 
tended discussion of these behaviors can be find in Ref.'"'^ . 

The Gaussian properties of the fluctuation amplitudes 
has also been investigated^^. Figure 17 reports the dis- 
tribution function of the resistance fluctuations p{R) ob- 
tained for a bias current / = 1.8 (A) (pre-breakdown 
region). The values of the parameters are the same used 
for the curves in Fig. 3. The parabola in Fig. 17 rep- 
resents the fit with a Gaussian distribution. The figure 
shows a strong enhancement of the probability for the 
positive fluctuations with respect to the Gaussian dis- 
tribution. This non-Gaussianity of the fluctuation am- 
plitudes increases at increasing current or at increasing 
temperature^^, when approaching the breakdown con- 
ditions. This emergence of a non-Gaussian behavior 
near the breakdown can be considered a relevant precur- 
sor of failure. Experiments performed in different ma- 
terials near the electrical breakdown confirm this non- 
Gaussianity of ffuctuations^^'^^'^^'^^. 

Finally, Fig. 18 displays the spectral densities of re- 
sistance fluctuations obtained for two increasing values 
of the bias current: / = 1.5 and / ~ 1.8 (A). The pa- 
rameters are the same as those in Fig. 17. The spec- 
tral densities have been calculated by Fourier transform- 
ing the corresponding correlation functions CsR{t). The 
spectra are found to be Lorentzian in agreement with 
experiments'*'^^'^^''*^. Moreover, the analysis of the corre- 
lation functions for different bias values and for different 
substrate temperatures shows that the correlation time 
is only weakly dependent on the bias while it is strongly 
dependent on the temperature^"^. 

IV. CONCLUSIONS 

I have presented a short review of a recently 
developed^"''*^'^^ approach which allows the study of the 
resistance noise over the full range of bias values, from 
the linear regime up to electrical breakdown. Resistance 
noise is described in terms of two competing processes 
taking place in RRNs. The two processes consist of the 
breaking and recovering of the elementary network re- 
sistors. The probabilities of the two processes are con- 
trolled by an electrical bias and by the external tem- 



perature. Monte Carlo simulations are performed to in- 
vestigate the RRN behavior as a function of the bias; 
steady state or electrical failure have been found. At 
the lowest biases, the two processes are practically ran- 
dom and an effective defect generation probability can 
be defined which controls the network behavior. In this 
Ohmic regime, a scaling relation has been found between 
the relative variance of resistance fluctuations and the 
average resistance^^. The value of the critical exponent 
is significantly higher than that associated with 1/f noise. 
The properties of the nonlinear regime, occurring when 
the bias overcomes a threshold value, are studied for dif- 
ferent values of the material dependent parameters and 
for different bias conditions (constant voltage or constant 
current)^"'^^. In general, two regions can be identified in 
the nonlinear regime: a moderate bias region, where both 
the average resistance and the relative variance of resis- 
tance fluctuations scale quadratically with the bias, and 
a pre-breakdown region where these quantities exhibit 
a superquadratic dependence. The quadratic behavior 
in the moderate bias region has been found independent 
of the values of the model parameters and of the bias 
conditions. A final conclusion about the universality of 
this behavior requires the investigation of the role of the 
network topology and of the electrical contact geometry. 
Remarkably, under constant voltage conditions, it has 
been found that the average resistance scales quadrati- 
cally over the full range of voltage values up to break- 
down. Moreover, it must be underlined the strong effects 
on the pre-breakdown properties of the temperature co- 
efficient of the resistance. Furthermore, non-Gaussian 
noise has been found in the pre-breakdown region under 
both bias conditions. The results obtained by simula- 
tions arc found to agree satisfactorily with electrical and 
noise measurements performed in composites^'^^'^^ and 
in semicontinuous metal films^^, and in the degradation 
processes of metallic lines due to electromigration^^'^^. 
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